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Abstract 

During the last fourteen years several exact results were ob- 
tained for the so-called boundary sine-Gordon model. In the 
case of a conformal bulk this 2D boundary quantum field the- 
ory describes the universal scaling behavior of the Caldeira- 
Leggett model of resistively shunted Josephson junctions. In 
this work, we use a recently developed Monte Carlo technique 
to test some of the analytical predictions. 
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1 Introduction 

For several decades, dissipation effects in quantum systems have attracted much interest, be- 
cause the coupling to the environment affects the properties of devices displaying macroscopic 
quantum coherence and may even induce phase transitions. These dynamical quantum phase 
transitions pose challenging theoretical problems (see e.g. book [1] and references therein). 
Among the simplest and most fundamental physical systems undergoing a dissipation driven 
transition is the resistively shunted Josephson junction, shown schematically in Fig.[IJ As 
long as the Josephson coupling energy Ej is much larger than the single electron charging 
energy Eq a superconducting state with a well defined phase difference across the junction 
is possible. In this case, the zero-temperature dynamics can be understood from the classical 
equation of motion for a single particle in a wash-board potential [2]: 

h 2 ■■ dUAS) ha ■ j i _ , hi , i , 

+ ^TT^ = -^> Ui{<t>) = -Ejcos{<t>)-— <f>, (1.1) 



8E C d(p 2tt t yT ' 2e 

where the dimensionless shunt conductance 

nh 



a = ^~2 R s 1 (1-2) 



plays the role of the damping strength. It readily follows from (II. ip that the zero-bias 
resistance, R = lim/^o ^4, vanishes (superconducting behavior). However, for so-called 
"ultra small" junctions characterized by Ej < Eq the phase difference in the superconductors 
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Figure 1: Circuit diagram for the shunted Josephson junction under an external current /. 
The junction is characterized by the Josephson energy Ej = ^- I c , where I c is the Josephson 

2 

critical current, and the single electron charging energy Eq = where C is the junction 
capacitance. It is shunted by the resistor Rs <C Rn, where Rn is the normal state tunneling 
resistance of the junction. 

exhibits large quantum fluctuations and the classical equation (11 .ip does not provide an 
adequate description of the system. Tunneling between local minima of the potential energy 
Ui(<p) (11. ip results in a finite voltage across the contact (insulating behavior) [3]. And because 
the effect of dissipation is to suppress the tunneling rate, the Superconductor-Insulator (S-I) 
transition in the ultra small junction is essentially controlled by the energy dissipation in 
the shunt. In the overdamped limit a — > oo, the phase difference <fi localizes in the vicinity 
of some minimum of Ui(4>) even for Eq 3> Ej and the device becomes superconducting [4]. 
On the other hand, for a = 0, tunneling events produce a voltage drop in sufficiently small 
Josephson junctions. There have been relatively little experimental studies of single, ultra 
small Josephson junctions. The a — ^ phase diagram for this important model system was 
measured in Refs. [5,6], and the experimental results are consistent with a zero-temperature 
S-I transition in the vicinity of the line a = 1 . 

The precise nature of the energy dissipation in real devices remains unclear. However, 
under certain conditions (see e.g. review [7]) the S-I transition is expected to be described 
by the model introduced by Caldeira and Leggett in their seminal work on macroscopic 
quantum tunneling [8]. At zero bias, the Matsubara action of that model is given by (here 
and below H = ks = 1) 

A = ll AT [i€ - E ""M 

While the first term in (11.31) readily follows from the classical equation, the last term in- 
troduces dissipation on the quantum level and results from integrating out a heat bath 
characterized by an Ohmic spectral density. As was argued by Schmid [9], the model (11.31) 
exhibits a dissipative phase transition at zero temperature, which occurs at a = 1 for suffi- 
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ciently small values of Jg-. Currently it is believed [10-12] that the boundary of the transition 
lies at a = 1 for any value of the ratio i.e., 

lim(A) = / ' ">1 , (L4 ) 

where R is the zero-bias resistance. 

A remarkable feature of the Caldeira-Leggett model is that it develops a universal scaling 
behavior characterized by a single energy scale 

E, oc E c (^A ^ , (1.5) 



provided 

a>l, (3E C ^oo, Ej/E c ^0. (1.6) 

In this regime the model can be studied by methods of Quantum Field Theory (QFT) and 
some physically relevant quantities can be calculated exactly [13-17]. 

Besides theoretical studies, numerical simulations of the model (II. 3p have been employed 
to investigate the location of the S-I phase boundary at intermediate Josephson coupling [18] 
and the scaling of the resistance with temperature and dissipation strength [19]. However, 
local update Monte Carlo (MC) schemes suffer from large autocorrelation times at low tem- 
perature, which prevent the simulation in the regime (11.61) . These limitations were overcome 
with the development of a cluster algorithm in Ref . [20] , which extends the powerful methods 
developed for spin systems by Swendsen, Wang and Wolff [21,22] to models described by 
non-local phase variables. The algorithm of Ref. [20] exploits the symmetry of the Joseph- 
son potential to build the clusters in such a way that no MC update is rejected and uses an 
efficient treatment of the long range couplings in the dissipation term [23]. Such non-local 
updates of the phase configurations dramatically reduce the autocorrelation time and allow 
an ergodic sampling of the discretized version of action (II. 3p . even at low temperatures and 
strong interactions. This recent progress in MC techniques provides an opportunity to ex- 
plore and test the relevance of the exact scaling results for the quantitative description of 
the Caldeira-Leggett model. The present study is an attempt in this direction. 

The main subject of our interest is the Matsubara two-point function in the Caldeira- 
Leggett model 

G,M = (e 1# We^(0)>, (1.7) 
at r = (3/2 and \q\ < |, in the experimentally relevant domain of parameters 

a>l, (3E C > 0.5xl0 3 , < 0.25 . (1.8) 

We begin by describing G g (f ) at zero temperature as a function of the dimensionless 
:ameters a, t 
Values (VEVs) 



parameters a, q and This limiting value is expressed in terms of the Vacuum Expectation 



limG a (f) = (e^>2 . (1.9) 

p — >oo 
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In order to clarify Eq. (11.91) one should note that in the domain a > 1 the global Z-symmetry, 

Z: 4>^4> + 27in (n = 0, ±1, ±2...) , (1.10) 

of the theory (11.31) is spontaneously broken (or equivalently, the associated quantum dissi- 
pative particle is localized at zero temperature). This implies among other things that the 
exponential operators gain nontrivial VEVs ( e 19 ^ )„, where by (•■•}„ we mean the expecta- 
tion value over the ground state in which is localized near Notice also that the first 
term in the expansion of ( e iq< ^ )o in q 2 gives the mean phase fluctuation at zero temperature. 
In Section[2]we combine an exact nonperturbative result from Ref. [17] with the conventional 
perturbative analyses to make quantitative predictions for the limiting value (II. 9p . Some 
technical details of the calculations are presented in the appendix. For = 0.125, 0.25, 
the analytical predictions for (jl.9p (q = 0.1, 0.25, 0.5) and for the mean phase fluctuation 
are found to be in good agreement with our MC results for the discretized action (I1.3P (see 
Figs.HandED. 

Unfortunately the problem of exact calculation of the scaling form of G q (^) has not 
yet been solved. In Section[3] we calculate the leading low-temperature correction to the 
scaling limit of G 9 (f ) in the case 1 < a < 2. It is proportional to (f3E*) 2 ~ 2a and cannot be 
actually used without paying attention to higher order corrections of the form oc (f3E*) 2k ( 1 ~ a ^ 
(k = 2, 3 . . .) (see Fig.H]). In the vicinity of the phase transition line a = 1 they can be 
accounted for within a Renormalization Group (RG) treatment. 

In SectionHJwe develop the proper renormalization scheme based on a certain choice of the 
running coupling constant. Namely, as it was pointed out in Ref. [24], the dimensionless ratio 
-j^- of zero-bias resistance and shunt resistance (or equivalently, the dimensionless mobility 
of the dissipative particle) possesses a universal scaling behavior, i.e., it is a function of the 
single scaling variable /3E* in the limit (II. 6p : 

/ R \ 

p = lim ( — ) (a>l). (1.11) 



l3E c -,oo \R S / 
0E* - fixed 

We shall argue that the scaling function p = p(f3E*) may serve as the running coupling (see 
Fig.[5]). This function was originally calculated for integer values of a = 2, 3, 4. . . within 
the Thermodynamic Bethe Ansatz approach in Ref. [15]. Later the closed expression for p 
was found for arbitrary a > 1 [16]. In Section H~2l we confirm the prediction from [16] by 
means of MC simulations (see Figs.[6]and[7]). Then, in Section H~3| we discuss the anomalous 
dimensions of exponential operators e iq< ^ and perform the RG re-summation of perturbative 
expansions for the two-point function G ? (|). The obtained results for G g (f) in the vicinity 
of the S-I phase transition line a = 1 have been confirmed by the MC simulations (see Fig. [8]). 

A recent numerical study [20] found that the S-I phase boundary in the Caldeira-Leggett 
model is a line of fixed points with continuously varying exponents. The universal behavior 
in this case is governed by a certain one-parameter family of boundary Conformal Field 
Theories (CFTs) whose modulus can be labeled by the limiting value 

p= lim ( — J : 0<p<l (<* = !)■ (!- 12 ) 



^c-» \R S 

Ej/Eq- fixed 



l It should be kept in mind that <f> in (|1 .3(1 is a non-compact variable — oo < (f> < oo, so that the different 
minima of the potential Uo(4>) — Ej cos(0) are distinguishable. 
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Nowadays this boundary CFT family is relatively well understood due to the works [25-28] 
which are mainly motivated by problems from String Theory (see e.g. [29]). In particular 
the exact value of the critical exponent r) q describing the universal behavior of the two-point 
correlation function ( 11. 71) . 

G q (r) oc \r\- 2r >" ( E c l < r < (3 ) , (1.13) 

can be extracted directly from results of Ref. [27] and reads explicitly as 

Vq = ^2 arcsin 2 (y/p sin(Trg)) (\q\ < \) . (1.14) 

In Section [5] we briefly recall some facts concerning the boundary CFTs associated with the 
critical line a = 1 and check the analytical expression for r) q by means of MC simulations. 
Our numerical results confirm the prediction (11 .14)) for q = 0.25 and 0.5 (see Fig.[HD. 

Finally, in Section[6], we briefly discuss the applicability of the scaling results in a wider 
domain of parameters than (II. 8p . As an illustration, we study the dependence p = p(§^) 
(11.121) along the critical line a = 1 for all values ^f- > (see Fig. fTTT) . 



2 Vacuum expectation values of exponential operators 



2.1 Small - expansion of (e^) 

It was realized a while ago [30] that the Caldeira-Leggett model in the scaling regime (11.61) 
can be thought of as an effective theory describing the dynamics of boundary degrees of 
freedom for the Boundary Sine-Gordon (BSG) model defined by the Euclidean action 

■Absg = j- / dx I dT [(d*V) 2 + (dr<p) 2 ] ~ 2p B / ^ dr cos(y/g<p B ) . (2.1) 

Here ip = ip(x, r) is a one component uncompactified Bose field, (x, r) are coordinates on 
the Euclidean half-cylinder, x < 0, t = t + f3, and 

<^(r)^(0,r) . (2.2) 

The interaction in (12. ID acts only on the boundary x = and is controlled by the two 
parameters g and ps- For < g < 1 the model is a CFT - a free Bose field with free 
boundary condition at x — - perturbed by the relevant boundary operator 2 cos( v / ^y? B ) of 
dimension g. We assume that this boundary operator is normalized such that the two-point 
correlation function takes the following asymptotic form 

<2cos( v ^^)(r) 2cos(^ i3 )(r , )> ^2|r-r'|- 29 as |r - r'| -> . (2.3) 
Under this normalization the parameter ps has the dimension [energy] 1-9 . 
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The original Callan-Thorlacius observation [30] is that upon fixing the boundary values 
of ifiB and integrating out the bulk part of the field (p, the action (12. ip reduces to the effective 
action (11.31) without the Coulomb term provided 

g = a~\ (2.4) 

and the variable is replaced by ^/gips- A superficial examination reveals that the model 
(II. 3p without the Coulomb term is perturbatively ill defined. It has to be equipped with an 
ultraviolet (UV) cutoff and needs renormalization. As long as the energy scale Eq remains 
finite but essentially exceeds all other energy scales, the Coulomb term in (11.31) just provides 
the desired regularization with the cutoff energy A oc ot.Ec- In what follows it will be 
convenient for us to specify the cutoff energy scale unambiguously, by 

A = aE c = 2.267 ... x aE c , (2.5) 

7T 

where 7^ = 0.5772 ... is Euler's constant. The dimensionless parameter 

E E 
e = — - = 0.220... x a" 1 — , (2.6) 
2A E c 

in turn will be treated as the bare coupling. Consistent removal of the UV divergence 
requires that the bare coupling be given a dependence on the cutoff scale such that the limit 
liniA^oo (A 1-9 e(A)) exists and coincides with the dimensionful parameter of the renormalized 
action (12.11) . One can always choose e(A) to satisfy the simple RG flow equation 

-A^ = (l-g)e 1 (2.7) 



V B = eA 1 - 9 . (2.8) 

The above relation immediately implies that Matsubara operators in the dissipative quan- 
tum mechanics (I1.3P can be thought of as bare boundary fields. Hence, they should admit 
asymptotic expansions in terms of the renormalized boundary fields from QFT (12. ip . In 
what follows we consider such an expansion for the exponential operators. 

The simplest and most fundamental boundary fields in QFT (12. ip are the exponential 
fields 

O q = e iq ^v* , 0\ = O g , (2.9) 

transforming irreducibly under the symmetry (11.101) : 

Z : O g -> e 2niqn O q . (2.10) 

Below, the field O q is always understood as a renormalized boundary field of definite scaling 
dimension 

d q = gq\ (2.11) 
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subjected to the UV renormalization condition 



T 



-2d a 



as 







(2.12) 



The general structure of the asymptotic expansion of the bare operator e 19 * readily follows 
from the standard RG arguments and from the global Z-symmetry (11.101) . (I2.10ln rl: 



e^(r) 



C^(q,e) A~ d ^ <V fc (r) + ... 



(2-13) 



fc=0,±l,±2. 



Here, the explicit dependence on the cutoff scale A is used to show the relative smallness 
of various terms and the dots in (12 . 1 3[) stand for the renormalized boundary fields of the 
form d T O q _k, d 2 O q ^k, (d T ip B ) 2 O q ^k ■ ■ ■ whose scaling dimensions are given by d q _k + N with 
N — 1, 2 ... . An obvious consequence of the <fi — > —(f) invariance of the action (12. ip is that 
the renormalization constants C^ k '(q,e) in (12.131) satisfy the relation 



C^ k \q,e)=C^ k \-q,e) 



(2.14) 



The renormalization constants do not carry any information on the infrared environment 
and therefore admit regular power series expansions in the bare coupling: 



C«(g,e) = el fc l J^cf (q) e 2 i 

3=0 



(2.15) 



The expansion coefficients are in principle computable within perturbation theory. In par- 
ticular, without tedious calculations, one finds 



a 



(0) 



(2.16) 



provided the cutoff scale is defined by Eq. (12.51) . The explicit expression for C^\q) 



Co (— <?) an d C{ (q) is somewhat cumbersome and relegated to the appendix (see Eqs. ( IA.8I) 
and flXT7p ). 

Contrary to the renormalization constants, VEVs of the boundary operators are impor- 
tant characteristics of the infrared environment. They do not possess a regular expansion 
in powers of e and, hence, can not be calculated perturbatively. In Ref. [17] the following 
explicit form for the VEVs of renormalized exponential fields (12.91) was proposed: 



(O q ) n = e 2 ^ n (p B )T=-° D(q) 



(2.17) 



where the function D(q) for |3?eg| < i is given by the integral 
D{q) = 



2% 



91 
1-9 



exp 



sinh 2 ( gqt) (e* 



2fl 
°° dt 



X 



1 + e 



t 



(2.18) 



+ e" 



-gt\ 



2 smh(gt) sinh(t) sinh((l — g)t) 



smh((l -g)t) 



+ e 



2 Although the Z-symmetry is spontaneously broken in the domain a > 1, the operator valued relations 
still possess this invariance. 

3 Here and below we use the symbol x to emphasize that the relation holds in a formal asymptotic sense 
without any reference to convergence issues. 
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and is denned through analytic continuation outside this domain. 

Combining now the exact result (I2.17P where /xg is given by (12. 8p with the asymptotic 
expansion (I2.13p . one obtains: 



jq<t> 



)n 



=0,±1,±2.. 



CW(q,e)D(q-k) 



min(l + 2g,2) ■ 

l + 0(e i- 9 ) 



(2.19) 



The remaining term 0(. . .) in (I2.19P requires some explanations. Namely, to transform the 
operator expansion ( 12.131) into the VEV expansion, small-e corrections to the BSG vacuum 
state should be taken into account. One can show that at the leading order this gives rise 

d q _ k +l +2g ' 



to the term oc e" 



1-9 



in (I2.19p . The correction oc e 



1-9 



in turn, is a contribution of 



the boundary field (d T (p B ) 2 O q -k of scaling dimension d q _k + 2 whose VEV is currently not 
known. 



2.2 Comparison with MC results 

For the numerical evaluation of ( e 1<?< ^ )o with \q\ < | we included only the contribution of the 
renormalized fields O q and O q ±i in the asymptotic expansion (12.191) . Furthermore, the per- 
turbative series for the renormalization constants C^ ±l \q, e) and C^(q, e) were truncated at 
the leading order (oc e) and at the second-order (oc e 2 ), respectively. Thus our approximation 
reads as follows: 



(e^) = [D(q) + Cf\q)D{q) e 2 + Cf (?) D(l - q) e 1+ ^ + 

CM(-q) D(l + q) e ^^ + 0( e mi ^^ 3+ ^))] < 1) . (2.20) 

Explicit formulas for the coefficients {q) and C^\q) can be found in the appendix. 
It deserves to be mentioned that for min(^-, j^-) < g < 1 and < q < | the last 
explicit term in expansion (I2.20p does not exceed the reminder 0(. . .) and therefore should 
be ignored]^ For the same reason the third term in (12.201) should be neglected if < g < 1 
and < q < |. Let us also note that the first nontrivial term of the expansion of (e 1<?< ^)o in 
a power series of q 2 is simply related to the mean phase fluctuation 



9(1-29) 



.,2 



:<**>„ = 1- \ <(0-0) 2 >o + O(g 4 ) , (2.21) 

where (f> = ( <fi )o = 0. Therefore Eqs. ( 12.201) and ( 12.2ip allow one to calculate ( (0 — 0) 2 ) at 

El 
e c - 



zero temperature for sufficiently small values of the ratio — 



4 Notice that it has poles in this domain of parameters. In order to clarify the general source of the 
singularities of the function D(q) (|2. 18|) . the definition of renormalized fields (|2.9| should be critically reex- 
amined. O q have been described as boundary fields with the definite scaling dimensions (|2.11|) and Z-charges 
(|2.10[) satisfying the renormalization condition (|2.12p . In general these requirements are sufficient to spec- 
ify Oq, except when certain resonance conditions are satisfied. We say that the renormalized field O^ 
has a n-th order resonance with the field O^ if the scaling dimensions of these fields satisfy the equation 
= +n(l- g), with some non-negative integer n. If this resonance condition is satisfied an obvious 
ambiguity, C*- 1 -* — > O^ + const (/is)™ O^, in defining the renormalized field O^ with the scaling dimen- 
sion typically results in the logarithmic scaling of O^. For example the simple pole of the function 
D(q) at q — just indicates that the field 0^_ has a first order resonance with O \__ x . 
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Figure 2: as a function of a 1 for q = 0.1, 0.25, 0.5 (from top to bottom). The MC 

data were obtained for = |, | with lattice spacing AtEc = \- The circles and diamonds 
correspond to /3J5c = 3750 and j3Ec = 7500, respectively. The dashed and solid lines are 
the theoretical predictions (I2.20p for the VEVs (e iq< ^)l for = | and ^ = ~, respectively. 
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The analytical prediction ( 12.20(1 has been tested by means of MC simulations performed 
on the discretized action 



X ' r ^ 1 a (4 - (t>jf 



i=0 



^ = " L SA rR. + ArEj C ° S(0l) + AW 2- -JT^Jh ' (2 - 22) 



8AtE c 



4N2 £j sin 2 (^' 



where iV is the number of time-slices of size At = (3/N . The algorithm of Ref. [20] combines 
local updates in the Fourier components 



N-l 



J2^ ink/N 4>n (2.23) 



n=0 



(required for ergodicity) with rejection-free cluster updates which can shift entire segments 
of phase variables from one minimum of the cosine-potential to another. 

The expectation value of (0 — 0) 2 can be obtained directly from the time domain variables 
<f> n . In order to measure the correlation functions (11.71) efficiently, it is useful to note that 



1 

G q (nAr) = - £ e~ 2 ™^ ( \^ k \" ) , (2.24) 



fc=0 
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Figure 3: The mean phase fluctuation as a function of a -1 . The MC data were obtained 
= |, | with lattice spacing AtE c = \. The circles and diamonds correspond to 
3750 and (3E C = 7500, respectively. The dashed and solid lines are theoretical 
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predictions for the VEV ( (<p — <p) 2 ) with J^- 



4 ana E C 



i 



respectively. 



where <& 9i fc is the k-th component of the Fourier transform of e 1<?< ^. By computing <& q ^ 
(in a time O(A^logA^) using fast Fourier transformation) and sampling |$ g ,fc| 2 in the MC 
simulation, one can therefore harness all the information contained in the phase configuration 
without considering each pair of phase variables. 

Results of the numerical work are depicted in Figs. [2] and [31 Similar results were found 
for = ^ and j|. No essential dependence of the simulation results on the discretization 
parameter AtE c was observed for AtE c ^0.5. 



3 Leading low-temperature correction for G q (^) 

We now focus on the considerable deviation between (e 19 ^)Q and the MC data for G 9 (|) in 
the domain g > 0.7 (see Fig.HJ). This discrepancy turns out to be a finite temperature effect 
and can be explored via infrared (IR) perturbation theory. 

To begin, let us briefly recall some well-known facts (see Refs. [14, 31, 32] for details) 
concerning the low-energy effective theory for the BSG model (12. ip . From the RG point 
of view this QFT describes a boundary flow between two fixed point. Each fixed point is 
controlled by the CFT of a free massless Bose field on the half-line x < subjected to a 
conformally invariant boundary condition at x = 0. Whereas the UV behavior corresponds 
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to the free (von Neumann) boundary condition, the IR fixed point is governed by the CFT 
with the fixed (Dirichlet) boundary condition. To be more precise, let H {n) be the space 
of states of the free massless Bose field on the half lane constrained by the fixed boundary 
condition tps = n| n {n = 0, ±1...). Then the whole space of states of the model ( 12.11) can 
be represented direct sum 



n 



(3.1) 



All terms in this sum are isomorphic to each other and the Hamiltonian associated with the 
IR fixed point, 

1 f° 

H* = — / dx (tt 2 + (d x ip) 2 ) where [ tt(x) , (fix 1 ) 1 = -2m 8(x - x) , (3.2) 
acts invariantly in TiS n \ 

H* : H (n) -> ft (n) . (3.3) 

The effective Hamiltonian governing the IR behavior can be introduced as a certain per- 
turbation by irrelevant boundary operators of the critical Hamiltonian ( 13. 2ft . To describe it 
explicitly we note that the set of primary boundary fields [33] of the infrared CFT includes 
{d T Cp) B = d T cp(0, t) and 



[e$*] 



B 



o 



m 



0, ±1, ±2... 



(3.4) 



where ip is the T-dua@ of (p. The boundary fields (13.41) intertwine the components in the 
direct sum ( 13.11) . 



and may be normalized by means of the condition 

(O m (T)oi(0)) =\T\- 2 -ir . 
With notation ( 13 Ah the effective Hamiltonian has the form 

H eff = H*-/i B (a i (0) + 6 t 1 (0)) + 



(3.5) 



(3.6) 



(3.7) 



Here dots stand for descendents of the identity operator of scaling dimensions 2, 4, 6 . . . 
acting invariantly in each linear subspace H,^. For our purposes the explicit form of the 
omitted terms is not essential, but we shall need an exact relation between the IR coupling 
/2b and the original parameter /xg in (12. ip [31,32]: 



2tt 



27T/i 



B 



Lr(i + ^)J 



r(i + a) 

2tt 



2vre 



r(i + a - 1 : 



A 



l-a 



(3.J 



°The T-dual of a free massless field is defined as usual, through the relations: d T <p ~ id x ip, d x (p — — i <9 r y. 
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Figure 4: G g (f ) as a function of a 1 for q = 0.1, 0.25, 0.5. The MC data were obtained for 
= \, (3E C = 1500 and AtE c = f . The solid lines are predictions (J22DJ for the VEVs 
(e'^Q. The dashed lines include the leading temperature correction (13.101) with /2g given 
by 



In the IR limit the exponential boundary fields (I2.12p should admit asymptotic expansions 
through the local boundary fields from the CFT associated with the IR fixed point, i.e., 



E 



^ D{q) I + D 1 (q) (d T Q>) B + 



(3.9) 



Here P n are projectors on the subspaces HS n \ I is the unit operator, and the dots stand 
for the contribution of the conformal descendents of I and (<9 T v?)# of scaling dimensions 
2, 3, 4 . . .. Unfortunately, the explicit form for the expansion coefficients in (13. 9p is currently 
not known. The sole exception is the function D(q) given by Eq. (12.181) . 

Using expansions (12.131) , (13.71) , (13.91) , it is now straightforward to calculate the leading 
temperature correction for G q (r) (11.71) . For r = | it is given by 



1 — 4 sin (nq) J(a) fi B 



f3 \2-2a 

2^) 



-min(2,4a— 4) 



(3.10) 



where the function J(a) is defined through analytic continuation of the integral 



J(a) 



du / dv 



4 sin 2 {^f) 



(3.11) 
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from the domain of convergence dtea < 1. The integral (13. lip can be expressed in terms of 
the generalized hypergeometric function 3 F 2 at unity: 



J (a) 



7r§r(i 



a) | 3^2(1,1 



a, 1; §,2 - a 1 1) 



2r(l - a) ' I- a 

3 -F 2 (l, 1 — a, 1 — a; I — a, 2 — a I l) 



(3.12) 



(1 -a)(l- 2a) 



It should be emphasized that the correction presented in Eq. (13.101) is the leading one for 
1 < a < 2 only. Indeed, the second term in expansion (I3.9P should give rise to a temperature 
correction oc (3~ 2 , dominating in the domain a > 2. In connection with this, one may note 
that the function J (01) (13.121) has a simple pole at a = 2. It is expected that the term oc f3~ 2 
cancels out this singularity and makes the expansion (13.101) finite at a = 2. Other important 
features of J (a) are that it vanishes at a = 1.3917 ... (g — 0.7185 . . .) and has a simple pole 
at a = 1. For this reason, one should expect that the leading low-temperature correction 
(13.101) can be applied literally in a narrow vicinity of a = 1.4. The MC results depicted in 
Fig. 0] are in agreement with this prediction. 



4 Gq(^) in the vicinity of the S-I phase transition 

One may argue that the low-temperature correction oc (/is) 2fc /3 2fc< - 1_a ) generated in the per- 
turbation expansion for the effective Hamiltonian (13.71) possesses a pole of order k at a — 1. 
For that reason the calculation of G g (f) for a close to one requires an effective summa- 
tion of all these singularities. Below we shall solve this problem using the RG improvement 
technique. 

4.1 Running coupling constant 

We begin with an appropriate choice of running coupling constant that provides a convenient 
global "coordinate" along the boundary flow in the BSG model including UV and IR fixed 
points. 

In general, the definition of the running coupling constant is ambiguous. As an example, 
the RG flow equation (12.71) suggests the following simple choice for a local coordinate in the 
vicinity of the UV fixed point: 

1 

Here E is the RG energy and E* oc /i^ -9 is the integration constant of the RG flow equation 
which sets the "physical scale" in the BSG model. At the same time the RG flow arrives at 
the IR fixed point along the irrelevant boundary field of scaling dimension a = g^ 1 . Hence 

/■*=(§) (4.2) 
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is another obvious choice for the running coupling. However, neither /uv nor /m can serve 
as a global coordinate for the RG flow, since they fail to describe the vicinity of the IR and 
UV fixed points, respectively. 

Notice that in this work we are dealing with the temperature dependence only, so that 
the RG energy E should be understood as the temperature and the dimensionless scaling 
variable 



K 



E* 
~E 



(3E* 



(4.3) 



is the inverse temperature measured in units of the "physical" scale. Also E*/p^~ 3 is a di- 
mensionless RG invariant which may be an arbitrary function of the RG invariant parameter 
g. It will be convenient for us to define E* unambiguously through the relation 



E. M 



1 


27T/i B 


1-9 A 


2vre 


2^ 


|r(i + $)J 


~ 2^ 


[r(i + s)J 



1 

l-g 



(4.4) 



It turns out that the scaling function p = p(n) (11. lip is a convenient choice for the running 
coupling for the boundary flow in the BSG model. As was mentioned in the introduction, 
this function is known exactly. It was originally calculated for integer values of a = 2, 3, 4 . . . 
within the Thermodynamic Bethe Ansatz approach in Ref. [15]. Later the closed expression 
for p was found for arbitrary a > 1 [16]. That result can be summarized as follows. Let \l/ 
be the solution of the ordinary differential equation 



-dl + s e 9y 



^ = 



such that 



(4.5) 



, y — > +oo 

const (y-y ), y -> -oo . 



(4.6) 



Then p can be expressed in terms of the function yo = yo(s) (14.61) : 

dyo 



p( K ) = 1 + g s 



ds 



provided 



2 2-2o 

s = g k y 



(4.7) 



(4-f 



With the exact result (14.71) one can prove (see [16] for details) that p admits the following 
high- and low-temperature expansions 



P (k) = 1-X>» (/uv) 2n 



^2p~n(u) {f 1R ] 



2n 



(4.9) 



n=l 



n=l 



Here the expansion in terms of the UV running coupling constant /uv = k 1 9 has a finite 
radius of convergence, while the power series in /jr = is an asymptotic expansion with 
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Figure 5: The RG flow for the BSG model in the (a, p) plane. The orientation of the RG 
trajectories (dashed lines) corresponds to a decrease of the RG energy E (temperature). The 
bold rays p = 0, 1 and the bold segment at a = 1 are lines of fixed points. 



zero radius of convergence. The coefficients in these series are functions of a satisfying a 
remarkable duality relation originally conjectured by Schmid [9]: 



p n (a) = p n {a 
The first coefficient p% explicitly reads [31] 



Pi (<* J 



,l-2as 



r 3 (i + a) 
r(| + a) 



(4.10) 



(4.11) 



As follows from Eqs. (I4.9p - fl4.lip . both /uv an d fm can be perturbatively expressed in terms 
of 1 — p and p, respectively, for any value of a. Hence the pair (a, p) may serve as global 
coordinates for the RG flow in the BSG model (see Fig.[5]). 

It is instructive to examine the RG flow equation in terms of the coordinates (a, p). 
Using Eqs. (I4.9P and (14.71) it is possible to show that 



E — 

dE 

E — 
dE 



. 

2(q-l)p(l-p) 
1 + {a - l)p 



1 - 



(a 



a 



U(a,p) 



(4.12) 



Here U(a,p) is a bounded functions for < a, p < 1 which admits the following convergent 
and asymptotic expansions 



U(p,a) = J2 U n 



la 



U(p,a) x J2u m (a) p r ' 



(4.13) 



m=l 



m=l 



Notice that the expansion coefficients satisfy the duality relation U m (a) = U m (a *). It is an 
immediate consequence of ( 14.121) that 



P = P + 0((a-l) 2 ) , 
6 Our numerical results suggest that < U(a, p) < min( 1, ( a "^s )• 



(4.14) 
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where p is the solution of the differential equation, 



dp 2(g-l)p(l-p) 
E dE~ l + (a-l)p ' (4 - 15) 

satisfying the asymptotic condition p — > pi (f3E*) 2 ~ 2a as E = f3~ l — > with pi given by 
(14.111) . One can easily integrate the differential equation (14.151) and show that p solves the 
equation 

p (i - p)-* = 2^ ^ Q^±4 mr*« . (4.i6) 

T(i + a) 

Despite the fact that the function p defined by Eq. ( 14.161) does not have a simple physical 
meaning, like the "physical" coupling p ( 11.111) . it may play the role of a running coupling in 
the same sense as p. In fact in the domain of our interest, i.e. 1 < a < 1.4, one can neglect 
the difference between p and p since the resulting numerical error turns out to be smaller 
than the error bars of the MC data. 

4.2 MC results for #- 



Numerically, we compute the zero-bias resistance by extrapolating the Fourier transform 
((f)(f)) UJk of (0(r)0(O) ) in Matsubara frequencies uik = — > 0: 



R a 
R~s = 2^ 



a 



J^o N 2 L 



^^72 1*1 (1^1 > „ . ( 4 - 17 ) 



where we used the notation (12.231) . It is therefore sufficient to sample the modules \(fik\ 2 of 
the lowest few Fourier coefficients. 

In Fig. [6] we compare the MC results for the ratio ^ as a function of (3Eq for a = 1.4 
and Ej/E c = \ against p (g2D and p fl4~T6j) with E* given by Eqs. (14.41) . (12.51) . Figure [7| 
shows the analytical prediction and numerical results for ^ as a function of a. Here, the 

MC data were obtained for (3E C = 1500, = | and AtE c = |, i.e., for the same values of 
parameters as in Fig.Hl 

Although the cluster algorithm allows us to study temperatures much lower than what 
was previously accessible (which means that the frequency points are densely spaced), dif- 
ferent extrapolation procedures can yield quite different results for R/Rs- A parabolic fit 
([2/0] Pade) to the lowest 9 Matsubara points gives the result shown by stars in Fig. [3 The 
downturn at the lowest frequencies is better captured by higher order Pade functions. We 
tried [2/1], [2/2] and [2/3] Pade approximations to the lowest 9 frequency points and found 
similar results for R/Rs- These estimates (shown as diamonds in Fig. [7]) agree reasonably 
well with the theoretical prediction for p and p. 



4.3 Anomalous dimensions 

Once we have chosen the global coordinates for the boundary RG flow, the scaling functions 
should be regarded as functions of (a,p). Here, our prime interest is the scaling limit of the 
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Figure 6: ^ as a function of (3E C - The MC data were obtained for a — 1.4, = | and 
AtE c = g. The solid and dashed lines represent p (14. T[) and p (I4.16p . respectively with E* 
given by Eqs. flO}, $U$, 



two-point correlator, 

lim G q (l) = G q (a,p) for \q\ < \ . (4.18) 

Let us first consider the case |g| < |. The renormalization theory predicts the general 
structure of the perturbative expansions of Q q . Namely, for |g| < | the exponential boundary 
operator is renormalized multiplicatively and, therefore, 

g^ t \a,p) = Z 2 q (a,p)W q (a,p) , (4.19) 

where the anomalous dimension, 

E dZr, 

and the factor W q possess (formal) power series expansions in p and 1 — p: 

oo oo 

W q (a,p) = J2 W * ( 1 "/ ? ) fc = zZ ibk P k ' 

fc=0 fc=0 

oo oo 



7 9 («,P) = ][> a - p)* = ^ • ( 4 - 21 > 



fc=0 fc=0 
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Figure 7: ^ as a function of a. The MC data were obtained for (3E C = 1500, §fe = \ and 
AtE c = | (the same as in Fig.HJ. The solid and dashed lines represent p fl4.T[) and p (I4.16p . 
respectively, with E* given by Eqs. 04.41) . (12.51) . (12.61) . 



It is important to bear in mind that both W q and 7 g are defined ambiguously, since the 
transformation W q — > W q F 2 , Z q — > Z q /F, i.e., 

preserves ^ pert ' ) . Here F is an arbitrary function in a rather broad class possessing a formal 
power series expansions at p = 1 and p = 0. This ambiguity does not affect the anomalous 
dimension at the critical lines p = 0, 1 and a = 1 (see Fig. [5]). At the line of UV fixed points 
p = 1, the operator O g becomes the primary conformal boundary field of scaling dimension 
d q (12.111) . On the other hand (as follows from Eq. (13. 9p ) O q becomes the unit operator in 
each subspace in the IR limit. Thus, we conclude that 

7 g (M) = ^, 7g(M) = 0. (4.23) 

Thanks to Refs. [27, 28] the value of the anomalous dimension at the S-I phase transition 
line a = 1 is known (see also Section[5] for some explanations): 

7<?(1,P) = -J arcsin 2 (y/p sin(Trg) ) . (4.24) 
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The inherent ambiguity (I4.22p in the definition of the anomalous dimension outside the 
critical region allows one to choose 7 g of the form 

7 q (a,p) = —j arcsin 2 (y/p sin(Trg) ) , (4.25) 



arr z 



which admits local expansions (14.211) and consistently fulfills the global restrictions ( 14.231) . 

All the above considerations suggest to resum the perturbative ( in / 2 R oc p 2 B /3 2 ~ 2a ) part 
of expansion ( 13.101) and, using Eqs. ( 14. 9ft . ( 14. IIP . express it in the form 



G<«"(f) 



(e"*)5 exp 



0:71" z 



[ — arcsin 2 (a/ p(n) sinfVg) ) 



w k p 



fc=l 



(4.26) 



where 



sin 2 (vrg) (2 2a ~ l V l 



a) 



0FT(l + a) 



J (a) + 



mat — 1, 



sin 2 (7rg) 



7T 



log4 + 0((a-l) 1 ) . 



(4.27) 



The major advantage of (I4.26[) is that all singularities at a = 1 are now absorbed by the 
exponential factor, i.e., the expansion coefficients Wk are regular functions of a in the vicinity 
of the S-I phase transition line. 

Similarly, to study G q (^) in the vicinity p = 1, it makes sense to write its high- 
temperature perturbative expansion in the form, 



1 + ^c fc (g) e' 

it=i 

oo 



exp 



2) Z^* d 



A 



arcsin 2 (\/ p{n) sin(7rg) ) 



x 



fe=i 



(4.28) 



Here the overall multiplicative constant should be chosen to satisfy the condition 



G< UV, (f) 



— 1d„ 



as (3 A > 1 > 







(4.29) 



which follows from the expansion (I2.13p . In particular, the first coefficient c\{q) reads 



Cl (q) = 2Cf\q) + 



47rgtan(7rg) \fnT{l + a~ l ) 



a- I 2Y(\ + a~ l ) 



(4.30) 



where C[°\q) is calculated in the appendix (see Eq. (|A. 17j) ) . Recall that the integration 
measure — in Eqs. ( I4.26[) . (14.281) is expressed through the running constant by means of the 
RG flow equations (14.121) . i.e., 



d/t 1 + (a — \)p 

~~k ~ ~2(a - l)p(l - p) 
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dp + 0((a- l) 1 ) 



(4.31) 



The values q = ±| require special consideration. In this case it is useful to start with 
the bare boundary fields 



cos 



(!) 



O- 



sin 



(I) 



(4.32) 



They carry the same charge with respect to the abelian group Z ( 11.101) . Nevertheless the 0± 
do not mix under RG transformations since they have different "C-parities" with respect to 
the flip <p — > —cf). Thus, the bare boundary fields 0± are renormalized multiplicatively with 
the corresponding anomalous dimensions 7±(a, p). Whereas these anomalous dimensions are 
the same at the critical line p — 1, 



7±(a, 1) 



1 

4a 



(4.33) 



they are different in the IR limit. Due to C-parity conservation one can predict that + 
and 0_ renormalize in the IR limit to the unit operator I and d T (p (in each subspace TC^ n ' 
(13. ip ). respectively, and hence 



7+(«,0) = 



7 _(a,0) = l . 



(4.34) 



An explicit form of j±(oc, p) at the critical line a = 1 follows again from results of Refs. [27, 
(see Section[5]): 



7+(l,p) = -5 arcsin 2 (v^) 



7T 



7 -(l,p) 



arcsin (y/p) 



7T 



Using Eqs. (14.331) . (I4.34p and (14.351) one can propose the following expressions for •y±: 

1 



7+(",p) 
7-(«»P) 



7i («, P) 



arcsin (-y/p) , 



(4.35) 



(4.36) 



1 arcsin (a/p) 

7T 



an' 



arcsin 



in 2 (Vp) • 



The correlator Gi(f) is expressed through the two-point correlation functions of 0± 



G UI) = G +(l)+ G -(l) > where = (O±(r)O±(0)) . 



(4.37) 



Notice that the anomalous dimensions 7 + and 7_ are well separated in the IR limit (see 
Eq. (I4.34p ) and the term G-(§) oc ((3E*)~ 2 in (14.371) is negligible in the low-temperature 
limit. For this reason, Eq. (I4.26P can be applied for q = ±= without any modification. At 
the same time the high-temperature perturbative expansion (14.281) is not applicable for \q\ 
sufficiently close to |. For q = ±| it should be replaced by (I4.37P where G± are understood 
as 



(I) 



exp 



fe=i 

oo 

i + r w w (i-p)i 



TT E; 

A 



7±(a,p(/c)) 



fc=i 



(4.38) 
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Figure 8: G q (%) as a function of a 1 for q = 0.1, 0.25, 0.5. The MC data were obtained for 
= \, (3E C = 1500 and AtE c = §. The solid lines are predictions (E2QD for the VEVs 
(e^)^ The dashed-dotted lines correspond to the RG-improved low-temperature expansion 
(14.261) . (I4.27p . The dotted lines represent the RG-improved high-temperature expansions 
(l4~28l) for q = 0.1, 0.25 and (O7ll - (IO0l for q = 0.5. All the terms in the perturbative series 
in 1 — p in ( 14.281) . except the first one, were neglected. 



It is straightforward to calculate the first nontrivial expansion coefficients in (14.381) : 

n 



' a— 1 ' 
2a , 



and 



2 



V^r(i- 

log4 - 

7T 



±2 -4 



_2V^r(i + a- 1 ) 



2a' 

- + - + 0((«-l) 1 ) 



7r a 



2a 2 
I) + ^ 



7T 



IT 



a — 1 



a 



2r(i + ari) J 



(4.39) 



(4.40) 



with C^\q) given by (1A.8I) . 

In Fig.[8]the MC data from Fig.HJare compared against the RG-improved low-temperature 
expansion (14.261) , (14.271) and the RG-improved high-temperature expansions ( 14.281) , (14.301) for 
q = 0.1, 0.25 and (I4.37I) - (I4.40I) for q = 0.5. Notice that all the terms in the perturbative 
series in 1 — p in (I4.28P were neglected except the first one. For the calculation of the integrals 
in (14T261) . (jQ8|) . (14381) we employed Eq. (I43il . Similar results were found for /3£ c = 375 
and 750. 
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5 Scaling dimensions of e iq( ^ at a = 1 



As follows from the RG flow equation (|2.7p . to achieve the scaling limit in the case g = 
a" 1 = 1 the bare coupling constant e (12. 6p should be kept fixed. Hence, the limiting value p 
(11.121) is a certain function of the dimensionless ratio It is discussed in the next section. 

Here we just mention that p is a monotonically decreasing function of (see Fig. fTT!) and, 
hence, the variable < p < 1 can be used to parametrize the critical line a = 1 to the same 
extent as J^-. Notice that the commonly used coupling for the BSG model at g = 1 is the 
one introduced in Ref. [26]. To relate p to the Calan-Klebanov-Ludwig-Maldacena coupling, 
gcKLMi one needs to repeat the calculation from the work [34] in the renormalization scheme 
of [26]. This yields 

\[P = COS (TCgcKLAl) ■ (5.1) 

The one-parameter family of boundary CFT associated with the line of fixed points a = 1 
possesses the global Z-invariance (11.101) . For this reason the corresponding boundary fields 
can be labeled by the "quasi-charge" 

-§<<?<§ (5-2) 

(analogous to the quasi-momentum of Bloch states in a periodic lattice), and the spectrum 
of their scaling dimensions splits into bands. Let us denote the primary boundary fields 
by 0(g, m ), where q is the quasi-charge from the first Brillouin zone (I5.2p and the integer 
m = 0, ±1, ±2... labels the different bands. According to the works [27,28] the scaling 
dimension of 0^ m ) is given by 



where 



lies within the range 



d(q, m ) = (\{q) +m) 2 , (5.3) 
X(q) = J arcsin [y/p sin(vrg)) (5.4) 

- i arcsin( v /p) < A < ~ arcsin( v /p) . (5.5) 

Thus, the allowed values of A(g) + m consist of bands of width - arcsin(y / p) centered at 
every integer, with gaps of width 1 — - arcsin( v /p). 

Similar to Eq. (12.131) . one should expect that the bare exponential operators for a = 1 
admit expansions of the form: 

e^(r) x C<m) (w) A ~ d( «' m) ( ^ m) {r) + ... , with q+\ = g + |(modZ) (5.6) 

m=0,±l,±2... 

and the dots stand for the conformal descendents of the scaling dimension d(q, m ) + iV (N = 
1, 2, . . .). The immediate consequence of (15. 6p is the expansion for the two-point correlator 



G 9 (r)x ^ An \~ 7 7 ) l Sm ~B ) ■ (5 ' 7) 

m=0,±l,±2... 
N=0, 1, 2... 



— 2dr/j ™i — 2N 
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Figure 9: The critical exponent Tjq clS 8b function of p (11.121) . The MC data were obtained 
for (3E C = 375, AtE c = §, j£ = 0.6, 1.0, 1.2, ... 2.0, and (3E C = 750, AtE c = f, j£ = 
0.2, 0.4, 0.8 (first, second and fourth point from the right). Errors on p are smaller than the 
symbol size. The lines show the predictions of Eq. (I1.14p . 



Here the dimensionless coefficients Affl are some functions of q and p. As long as \q\ < | 
and (3 A 3> 1, the first term with m = N = dominates in this asymptotic series: 



G,(r 



( T ) „ 4»> 



\ TT 



sin 



7TT 







l?l < I) 



(5.f 



where % = c?( 9)0 ) fll . 14[) . Notice that, as follows from Eqs. (l5.3p . (15 .4p . c?(i,o) = ^(±,-1) f° r 
p = 1. Hence, near the boundary q — | of the first Brillouin zone (15.21) and for p close to 
one, it makes sense to keep also the first subleading term with m = —1, N = in (15 ,7p . In 
particular, 



Gijr) 



sin 



7TT 







-2,, 



+ 4 



sin 



7TT 







-2r/ 



(5.9) 



where 



(!-v^) 2 - 



(5.10) 



One may note that the limit a — > 1 brings the UV expansions (14.281) and (14.371) . (14.381) into 
the forms (15.81) and (15. 9h . respectively, taken at r = /3/2. 
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Figure 10: Coefficient v ( 15. lip as a function of p ( 11.12)) . The points were obtained by fitting 
the MC data for f3E c = 375, 750, 1000, 1250, 1500 with AtE c = § and (from right to left) 
g = 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4, 1.6, 1.8, 2.2, 2.6, 3. Errors on p are smaller than the symbol 
size. The solid line is v = p. 

We have performed the numerical check of Eqs. (15. 8 p and (15.91) . (I5.10p . Namely, we used 
these general forms to fit the function G q (r) obtained from the MC simulations. The con- 
stants Aq ^ and r] q were treated as fitting parameters. In Fig. [9] the results of these 
fits for i] q are compared against the analytical prediction (11.141) for q = 0.25, 0.5 and several 
values of p. 

We have also numerically studied the mean phase fluctuation for a = 1. As it was 
observed in Ref. [20], the phase fluctuations grow proportional to the logarithm of the 
inverse temperature, 



In Fig. [TD] our numerical results for z/, obtained by fitting the form (15. lip to the MC data for 
( {(f) — <j)) 2 ), are plotted against p. In all likelihood, 



( (0 - 0) 2 ) -> 2u log(/3£ c ) + const 



{fiEa - oo) . 



(5.11) 



v = p 



1 d 2 d q . 

2 ~£hf 



q=0 



(5.12) 
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6 Nonperturbative energy scale 



Scaling functions like p (11. lip and Q q (14.181) are defined as a particular limit of their cor- 
responding "bare" counterparts from the Caldeira-Leggett model. However, as has been 
demonstrated above, they can be successfully employed for numerical approximations of the 
bare functions for nonzero e and finite values of (3A provided the energy scale E* is under- 
stood as in (14. 4p . Of course, the bare functions differ from their scaling limits by corrections 
vanishing as e — > and the quality of the "scaling" approximation diminishes as the bare 
coupling e is increased. It turns out that some of these corrections can be captured, and 
the overall quality of the scaling approximation considerably improved, by substituting E* 
with some nonperturbatively defined "physical" scale for the Caldeira-Leggett model. Recall 
that E* has been introduced as a first integral of the RG flow equation (12. 7p . i.e., it is a 
perturbatively defined RG-invariant energy scale in the BSG model. 

( CL*) 

To introduce the physical scale in the Caldeira-Leggett model, El , we take advantage 
of the result of Ref. [35]. Namely, Weiss and Grabert found that at low temperatures 



R s ~ 2T{\ + a) V 7T J \ujj 
where 

VZEjEc (6.2) 



and A is the so-called tunneling amplitude (see also book [1]). It is important that the valid- 
ity of the leading low-temperature asymptotic ( 16.11) does not actually require the smallness 
of the ratio and therefore, the nonperturbative energy scale can be identified as 



E (CL) = 

2tt 



7rA 



(6.3) 



//-IT \ 

It is easy to see that Eq. ( 16. ip , when expressed in terms of El , has the form of the leading 
low-temperature correction to the scaling function p (14.91) . (14.111) provided E* is replaced by 

//-IT \ (C^X \ 

El . One should expect that the scaling function p (14.71) with n replaced by (3 El well 
approximates -j^- even for jj£ > 1, provided the temperature is not too high compared to 
the cutoff scale A cx aEc- By virtue of Eqs. (12.51) . (12.61) and (14. 4p . one has 



E c 



2T(1 + a 



vr 2 e" 7B 



4r(a~ 



[zF(z,a)] aa -\ (6.4) 



where F is some function defined for any 



z=|i>0, (6.5) 

such that F(0, a) = 1. 

The function F = F(z,a) in Eq. (16.41) can be studied using conventional perturbation 
theory in the bare coupling. In particular, at the second order in perturbation theory we 
find 

F(z, a) = 1 + h(a) z 2 + O ( z min(4 ^ } ) , (6.6) 
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where 



/i(a) = §Cf(l) (^)\ (6-7) 
with c[ (1) the coefficient (g) ( 1A.17I) taken at g = 1. Notice that has a simple 

a + 2 

pole at a = 4 which is expected to be canceled by the nonperturbative correction oc Z"- 1 , 
and, therefore, the term presented in Eq. (16 .6p is a leading small- z correction in the domain 
1 < a < 4. Numerical values for fi(a) for 1/1.2 < a < 2 are presented in Tabled] 



9 


a 


fi 


0.50 


2.00 


0.11513 


0.60 


1.67 


0.13243 


0.70 


1.43 


0.15176 


0.80 


1.25 


0.17177 


0.90 


1.11 


0.19203 


1.00 


1.00 


0.21234 


1.10 


0.91 


0.23260 


1.20 


0.83 


0.25278 



Table 1: Coefficient f\ (16. 7ft as a function of a = g 1 . 



In the strong coupling limit, — > oo, the energy scale Ei Ch ^ can be explored by means of 
the saddle-point approximation and the original instanton calculation from Ref. [9] suggests 
the following semiclassical form for the tunneling amplitude A in ( 16. 31) : 



A 

2oo, 



iz 4 e 



0(z-*) 



(6.8) 



To illustrate the results from this section, let us consider the -Ej/^-dependence of p 
at a = 1. Due to the relation (I4.14p . one can apply Eq. (14. 161) with E* replaced by Ei CL ^ 
( 16. 3p . (16 .4p as a — > 1. Then, for a = 1 one finds: 



1 



where X 
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zF(z,l) . 



(6.9) 



1 + X 2 ttA 4e^ 

We may now apply the weak and strong coupling expansions (16.61) . (16.81) to the case a = 1: 



P 



(weak) 



1 + X 2 



X 



TV 



[l + f 1 z 2 + f 2 z i + ...] ( 2 = g«l), (6.10) 



and 



P 



(strong) 



57T K P e 



1 + az- 1 z +0(z~ 1 ) 



( 



The numerical value of the coefficient f\ in (16.101) is presented in Table[TJ 



a=l 



0.21234. . . 
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(6.11) 



(6.12) 



E.j 
E c 


P 


n 
r 


{ weak) 


(strong) 

r 


0.2 





918 ± 0.004 





927b 




A A 

0.4 





753 ± 0.005 


A 




753 




A r> 

0.6 





565 ± 0.004 


A 




C CT O 

558 




A o 

(J. 8 





392 ± 0.004 


A 




OA < ' 

396 




1 A 

1.0 





275 ± 0.004 





278 




1.2 





186 ± 0.003 





OA 

20 




1 A 

1.4 





126 ± 0.002 





15 


0.15 


1.0 





086 ± 0.003 





"1 O 

12 


A AAO 

0.098 


1.8 





0622 ± 0.0015 






0.066 


o a 

2.0 





A A OA 1 A AAO 

0420 ± 0.002 






A A A P 

0.046 


2.2 





0315 ± 0.0008 






0.032 


2.4 





0226 ± 0.0008 






0.0228 


2.6 





0168 ± 0.0005 






0.0164 


3.0 





0090 ± 0.0004 






0.0088 


4.0 





0022± 0.0002 






0.0022 



Table 2: p fll . 1 2j) as a function of for a = 1. The MC data were obtained for (3Ec = 375 
and AtE c = f, ^. The expansions p (weak) and p( stron s) are given by flBTTUD and fl6TTT| . 
respectively, with the numerical coefficients chosen as f± = 0.21234..., fi w —0.05 and 
a w 0.5. 

Coefficients f'2 (I6.10p and a (16. lip are currently not known, and we have estimated their 
values using the available MC data (see Tableland Fig. fTT]) : 

f 2 w -0.05 , a w 0.5 . (6.13) 



7 Conclusion 

The resistively shunted Josephson junction, represented by the Caldeira-Leggett action with 
non-compact phase variable, is one of the simplest and most fundamental models exhibiting a 
dissipation driven quantum phase transition. It has been studied theoretically for more than 
twenty years. Most recently, there has been considerable interest in phenomena attributed 
to the interplay between spatial couplings and dissipation in systems involving several junc- 
tions [36,37] and Josephson junction arrays [38-40]. Complicated phase diagrams and lines 
of fixed points have been predicted for the latter models using strong- and weak-coupling 
approximations. 

The purpose of the present paper was to highlight some analytical results for the single 
junction model, which - at least in the solid state community - have not been widely ap- 
preciated. We took advantage of the theoretical machinery developed in the QFT/String 
Theory community to discuss analytical expressions characterizing the behavior of the junc- 
tion at criticality and in the superconducting phase. In particular, we studied the Matsubara 
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a. 




Figure 11: p fl!.12[) as a function of -f^- for a = 1. The MC data were obtained for 
(3E C = 375 and AtE c — |, ^. The solid and dashed lines correspond to weak coupling 
(16.101) and strong coupling (jfi.lip expansions, respectively, with the numerical coefficients 
chosen as f x = 0.21234 . . ., f 2 « -0.05 and a « 0.5. 

correlator (e 19 *(r) e - " 39 ^) ) as a function of a, g, Ej/E c and temperature. On the S-I phase 
boundary, this correlator decays with a critical exponent, which is a function of the dimen- 
sionless zero-bias resistance p. We also discussed an implicit expression for p as a function 
of a and temperature. The second quantity of interest to us was the mean phase fluctuation 
((0 — 0) 2 ). At the critical line a = 1, the fluctuations grow proportional to the logarithm of 
the inverse temperature, ((0 — 0) 2 ) — > 2u\og{j3E c ) + const as (3 — > oo, and we have proposed 
that v = p. 

The accuracy of the analytical formulas and the range of applicability of the perturbative 
calculations have been tested by means of MC simulations. These numerical checks were 
made possible by the recent development of an efficient cluster sampling scheme [20]. Hence, 
the remarkable agreement between analytical and numerical results not only verifies the 
accuracy of the formulas discussed in this paper, but also demonstrates the performance of 
the cluster algorithm. 
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A Appendix 

In the case of the Caldeira-Leggett model, the coefficients in the perturbative series ( 12.151) 
do not possess a simple analytic form. Here, we consider integral representations for cf > \q) 
and Cq ±1 ' ) (q) which are sufficient for their numerical evaluations. 

The perturbative propagator for the action ( 11. 3ft with Ej = explicitly reads 



/ , / \ , /„x v 47rg cos(co> n r) . . _ . 

( T pert = -f Yl ■ ( A ' X ) 

Here uj n = =^ are Matsubara frequencies and notations (12 ,4p . (12. 5p are applied. As the 
temperature goes to zero, one has 

(0(r)0(O)) pcrt = 2^1og(|f)+log0(Ar) + O(^) , (A.2) 

with 



<&(v) = exp 



sin 2 (^ 
Ag / dv 2 



o 



(A.3) 



v + v 2 e^ E 

Notice that the zero-temperature Green function ( 1A.3I) satisfies the conditions 0(0) = 1 and 
(5(u) = \ v \- 2 9 (i + o( v - 2 )) as |u|->oo. (A.4) 

For $te(qg) > | the perturbative coefficient Cq (g) is given by 

/>oo 

C^g) = dv <S q {v) . (A.5) 



The coefficient Cf^ (q) also admits a simple integral representation 

/oo poo 
dv / du[(5(?;-u)<5 9 (?;)<5- 9 (u) -0( V - M )] (A.6) 
■oo J — oo 

within the domain of convergence 

5 > 1 , -q < s Re (q) < q with q = 1 - ± . (A.7) 

Both integrals (IA.5I) . (IA.6I) diverge in the domain of our interest, i.e. for < g < 1 and 
\q\ < |. Such divergences are of a similar nature as those encountered in calculations of 
the structure functions within conformal perturbation theory [42]. Namely, they indicate 
the presence of nonperturbative terms in the e-expansion (I2.19p . i.e., this is a price we have 
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to pay for the analyticity of the renormalization constants (I2.15p . To assign an appropriate 
meaning to expressions (1A.5|) and (1A.6|) in the domain of interest one should regularize the 



integrals somehow and then extract their finite parts. Of course, the finite parts are not 
fixed unambiguously by any general principle. In a sense, this implies an ambiguity in the 
choice of the bare coupling. As is usual in QFT, in the absence of logarithmic divergences 
(which is the case) the analytical regularization (an analytic continuation in field dimension) 
proves to be the most convenient for calculations. In the case under consideration it implies 
a certain renormalization scheme based on the simplest admissible RG flow equation (12.71) . 

For the above reasons we introduce (q) and (q) through an analytic continuation 
in g (the scaling dimension of the perturbation in (I2.ip ) of the integrals (1A.5I) . flA.7j) from 
their domains of convergence. In the case of C^lq) the analytic continuation can be easily 
performed by rewriting (1A.5I) in the form 

poo 



where 



and 



C™(5) = / d« [<S'M - «'(«) ] + <T(9) , (A. 



c* iq) = (A . 9) 



e(v) = {l + v 2 )- g . (A.10) 



Now the integral in (1A.8I) converges for 9fte(gg) > — ~ and defines an analytic function of 
the variable gq. This implies in particular that Cq^(1 + g), considered as a function of the 
complex variable q for < g < 1 has a simple pole at q = —go = ^ ~ 1 with the residue: 

C$ 1) (l + ?)--^- as q^-qo. (A.11) 

q + qo 

Notice that the coefficient C$ ^ (q) can be calculated using ( 1A.8I) and the relation Cq ^ (q) = 

An analytic continuation of the integral (1A.6I) is slightly more subtle. Before proceeding 
with its explicit description let us make an obvious remark concerning the function flA.lOj) . 
It is easy to see that log<S(Ar) can be thought of as a perturbative, zero-temperature prop- 
agator for the model defined by the Matsubara action 

i 00 i i /■! 

A = ^T1 S 1^1 e ~ \<t>{n)\ 2 -Ej I drcos(0), (A.12) 

n=-oo 

I 

where <f)(n) = f 2 ? dr0(r) e luJnT . Despite the lack of a physical interpretation, the action is a 

suitable UV regularization of the boundary effective action for the QFT (12.11) . Furthermore 
it has pleasant advantages compared to the Caldeira-Leggett action ( 11. 31) . Namely, the first 
perturbative coefficients in the expansions of renormalization constants for ( 1A.12I) (which 
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will be denoted as Cj^\q)) admit relatively simple analytic forms. The coefficient C^lq) is 
given by COD, while C[°\q) can be expressed in terms of the generalized hypergeometric 
function 4F3: 

2 3 -M%!r(-i + (i- g ) ff ) 



3%) 



X 



cos(7rg#)(l -g)T(qg) 
4F^(-l-l + g,-qg,-l + (l-q)g; ± - qg , -± + ^ g , 1=2 g | i) 
2 i-(i+ 9 ) 97r | r( _i + (l + q ) g ) 



x 



where 



cos(7rg#)r(-g#) 
4F^{ 1 1 ,g,qg,- 1 1 + (l + q)g; ±* g, § + g<?, | + p | 1 ) 

1 

4 F 3 reg (a,6, c,d;ej,g\z) 



(A.13) 



r(e)r(/)rfo) 



(A.14) 



. ( afecrf 2 a(a + 1)6(6+ l)c(c+ 1) 
+ H + e(e + l)f(f + l)g(g+l) 2! + 



To write down an efficient integral representation for C[°\q) in the Caldeira-Leggett 
model we first observe that C[°\q) has three simple poles located at the boundary of the 
domain of convergence (1A.7j) : 

27rgtan(7rgg) 



as g — > 1 



and 



C [°\q)^C^(lTq)C^(±q) as ±q 



,(i) 



1 - 
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(A.15) 



(A.16) 



where Cq (g) is defined by Eq. fl A.8|) and satisfies flA.111) . This observation implies that for 
numerical evaluation of C[°\q) = Ci(—q) in the domain < q < 1 and g> > it makes 
sense to rewrite (1A.6I) in the form, 



where 



5C[°\q) = CF(-q)C?>(l + q)- C^(-q)C^(l + q) + 



(1), 



,(1), 



dv 



(A.17) 



(A.18) 



a 



(1) 



g)(5» - C^(-q)& 1+(1 (v) + ^(-g)© 1+(? (t;) - 0^(1 - q)&{v) + 

/ du [ ®{v - u) <5 q {v)<&- q {u) - <S{v -u)- <S q {v)&- q (u) -(0^0)]}. 
./ —00 

It should be noted that for the numerical results depicted in Figs. [2] and [3] it is sufficient to 
approximate C[°\q) by 

27rgtan(7rgg) 



+ Ci%)C£>{l -q) + (%>{-q)C?>(l + q) . 



1-5 



(A.19) 



The numerical error due to this approximation is considerably smaller than error bars of our 
MC data. 
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